Impact of cosmic inhomogeneities on SNe observations 
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Abstract. We study the impact of cosmic inhomogeneities on the interpretation of SNe observations. We build an inhomoge- 
neous universe model that can confront supernova data and yet is reasonably well compatible with the Copernican Principle. 
Our model combines a relatively small local void, that gives apparent acceleration at low redshifts, with a meatball model that 
gives sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different effects of voids on 
the data, allow the model to mimic the concordance model. 
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INTRODUCTION 

Until the nature of the dark energy is completely understood, the "safe" consequence of the success of the concordance 
model is that the flat isotropic and homogeneous ACDM model is a good phenomenological fit to the real inhomoge- 
neous universe and not that the concordance model is the actual model of the universe. It is therefore useful to look 
for alternative models that fit the data. Here we will discuss the possibility that the late time evidence for dark energy 
can be explained by the late time formation of nonlinear large-scale inhomogeneities, so that the so-called coincidence 
problem becomes a hint of new physics. 

In the following, we will assume that the spacetime of the inhomogeneous universe is accurately described by small 
perturbations around the Friedmann-Lemaitre-Robertson- Walker (FLRW) solution whose energy content and spatial 
curvature are defined as Hubble-volume spatial averages over the inhomogeneous universe. Following Ref. [3] we will 
call this Hubble-volume average the Global Background Solution (GBS). The cosmological model obtained through 
observations, on the other hand, will be called the Phenomenological Background Solution (PBS). The idea is to 
associate the concordance model with the PBS, while the actual GBS is the Einstein-de Sitter (EdS) model. 

SETUP 

For consistency with the CMB spectrum and the age of the universe, the current background expansion rate of the EdS 
model is taken to be H„ = 100 h m km s _1 Mpc -1 with hoo = 0.5. The local expansion rate Hq will be higher, about 
h ~ 0.58 in the current model, which is within 2-0 of the HST key project [4] value of 0.72 ± 0.08 JH]. 

Our setup is made of two elements. First, we will model the overall universe by a meatball model consisting of two 
families of randomly placed halos. The family A describes very large, low density contrast structures and the family 
B models large clusters of galaxies. The parameters specifying the model are (for each family) the average comoving 
distance between meatballs d c , the proper radius of the meatball R p and the mass of the meatball M. The numerical 
values we used are given in Table Q] 

The distance d c between meatballs is connected to their comoving number density by n c = K J 6 7 3 . We assume 
the meatballs to have virialized at z = 1.6 and therefore their proper radii R p are constants. For the family B we 
used the singular isothermal spheres (SIS) density profile and defined the radius to have a density contrast of 200 at 
virialization. For the family A the gaussian profile was used and the radius was chosen to qualitatively describe the 
largest structures seen in simulations of structure formation. 



1 Contribution to the proceedings of the Invisible Universe International Conference, Paris, France, June 29-July 3, 2009. Based on |1, 2], where 
more details and references can be found. 



TABLE 1. Parameters of the meatball model (qq = !)■ 



Quantity 


Family A 


Family B 


(k 


100/T 1 Mpc 


10 hT x Mpc 


R P 


\0h- 1 Mpc 


580/1" 1 kpc 


M 


6.1- 10 17 /r _1 Af© 


6.1 • 10 14 A _1 M e 


density profile 


Gaussian 


SIS 



The distribution of family A meatballs leaves underdense regions of order lOO/i Mpc which are filled by the 
family B meatballs. Because the mass is equally subdivided between the two families, these underdense regions have 
on average a density contrast of 8 « —0.5. The second element of our setup consists of placing the observer in one 
of the large underdense regions and modelling the metric of this particular void around the observer more accurately 
with a Lemaitre-Tolman-Bondi (LTB) bubble matched to the EdS background metric. See Fig.Q]for a sketch. 




FIGURE 1. The larger disks represent the meatballs of the family A, while the smaller ones represent the meatballs of the family 
B. The concentric circles mark the overdense shells in the LTB bubble and are roughly 100/z~'Mpc apart. 



LENSING 

In a late-time universe dominated by voids the homogeneity is recovered only on scales larger than the largest inho- 
mogeneity scale, which in our model is L/,„„, ~ lOO/i^'Mpc. However, type la supernovae (SNe) probe angular scales 
LsNe ^ ^hom an d it is therefore not clear that the physics inferred from SNe observations can be directly associated 
with the smoothed-out GBS model, especially with small data samples. Epigrammatically, the commutativity between 
averaging and measuring is not guaranteed. Note that this idea belongs to weak-backreaction studies and is differ- 
ent from the possible non-commutativity between averaging and dynamics which is the kernel of strong-backreaction 
studies (see for example Ref. [0]): 

[(•••),EoM] ^0 => ABS^GBS 
[(•••), Obs] ^0 =S> PBS + GBS 

See JH for more details about a proposed framework for the backreaction proposal. 

Cumulative gravitational lensing is a possible source for the non-commutativity between averaging and measuring 
and to study this effect we have focused on a meatball model. The latter indeed incorporates quantitatively the crucial 
feature that photons can travel through voids missing the localised overdensities. This feature, instead, is not present 
in swiss-cheese models where the bubble boundaries are designed to have compensating overdensities. Such models 
have indeed been shown to have on average little lensing effects |0] • 



FIGURE 2. Shown is a co-moving segment of a photon geodesic between an observer at O and a source at S. The shaded disks 
represent meatballs of radius R. 



In the weak-lensing approximation the lens convergence K is given by 



K(z)= drG{r,r s {z))8{r,t(r)), (1) 



where 5(r,f) is the density contrast and 

3Hj r(r s -r) 1 

The functions a(t), t(r) and r(z) correspond to the EdS background model, r s is the comoving position of the source at 
redshift z and the integral is evaluated along the unperturbed light path. We will denote by Kg the convergence of the 
"empty beam" which is given by Eq. (fTJ with 8 = — 1 . The shift in the distance modulus, neglecting the second-order 
contribution of the shear, is 

Am(z)=5lDg 10 (l-ir(z)). (3) 

Am is proportional to the total matter column density along the photon path. For a lower-than-GBS column density the 
light is demagnified, while in the opposite case it is magnified. 

In Ref. |01 we derived a fast and easy way to obtain the convergence probability distribution function (PDF) for 
meatball models. First define a normalized profile <p(x) = p(x)/p, where p is the density profile of a meatball and p 
is the EdS density. Then the following quantity characterizes the contribution from one meatball, hit with an impact 
parameter b, to the convergence (see Fig.|2]for an illustration of the calculation): 

, v r R ® 2xdx , 
T(b,t)= ^_^ y(jt,0. (4) 

Jb Vx 2 - b 2 

Then divide the comoving distance r s to the source and the radius R of the meatball into Ns and bins of widths 
Ar, and Ab m , such that Ab m <C R and R <C Ar, -C r s . The convergence due to a meatball placed within the bin (i,m) is 
just 

Klim = G{n,r s ) T(b,„ ,t t ). (5) 
The PDF is then generated by the following functional 

K ({kim}) = L L m ' <® 

,= lm=l N 



where k[NoAN]im is a Poisson random variable of parameter NoANi,„. No[z) is the number of observations at redshift 
z and ANj„, = n c AV, m is the expected number of meatballs in the volume AV, m = 2nb m Ab m Arj. For a Poisson variable 
of parameter A , the mean is A and so the expected convergence in Eq. © is zero, consistent with photon conservation 
in weak lensing. 

The PDF is generated by creating a large sample of configurations {£,-,„} drawn from the Poisson distribution, and 
computing the appropriate convergences through Eq. (O. This is a task that a laptop fullfils in less than a second as 
opposed to expensive ray tracing techniques. The weak lensing approximation of Eq. (O describes the mode of the 
PDF very well always and, in case of non-point-like meatballs, also the tails of the PDF with < 5% of error. See Fig. [3] 
for an example relative to the present model. 
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FIGURE 3. Lensing PDF in magnitudes for the meatball model of this paper at redshift z = 1 and for a set of 1 and 10 SNe. 
Increasing the number of measurements reduces the skewness of the PDF, which, however, is still highly non gaussian after 10 
measurements. 

Eq. © displays explicitly the effect of the size of the data sample: even if K might have a skewed PDF and a nonzero 
mode for No = 1, the gaussianity is recovered for large No and the set of observations becomes unbiased. To define 
No{z) we will use the Union Catalogue of Ref. @l, which consists of 307 SNe spread between < z < 1 .6. We have 
binned the SNe with a bin width of Az = 0.1 and the result in plotted in the top panel of Fig. |4] SNe observations 
within the same bin are then be treated as repetition of the same observation. 

An analytical approximation for the PDF was found in Ref. (2l by resumming the individual Poisson varaibles. The 
result is simply 

where k is again a Poisson random variable and 

N = n c -r s {z)E-nR 2 Q\. (8) 

Physically N corresponds to the expected number of collisions with metaballs in a tube of a radius RQy and a length 
r s E connecting the observer to the source. The constants E ~ 2/3 and 2<p are corrections for the effective comoving 
distance and the effective meatball radius. For the gaussian meatballs of the first family ~ 0.5 and for the SIS 
meatballs of the second family ~ 0.25. The co-moving size R is computed at half distance between observer and 
source. Eq. (0 reproduces the results of Eq. (O quite well for a wide range of parameters, and so Eq. ([8]l provides an 
easy way to estimate the magnification bias. If meatballs fill the universe smoothly (large n c ) or their density profile is 
wide and smooth (large R) then N is large leaving a negligible magnification bias. Moreover, even for a small N a large 
No removes the non-Gaussianity from the effective PDF. However, if the product NoN is small, the lens convergence 
approaches the empty beam limit as the light rays "scan" mostly voids. This gives rise to a low effective column density 
and therefore to an average demagnification. For a thorough analysis of the previous results we refer to Ref. 0]. 
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FIGURE 4. Top: the 307 SNe of the Union Catalogue of Ref. [ 8] binned with Az = 0. 1. Bottom: shift in the distance modulus 
Am for Nq taken from the top panel and for the meatball model of this paper. 



We have used a modification of the turboGL package |0] to evaluate the magnification bias for our model and 
for the Union Catalogue. Results are shown in the bottom panel of Fig. 2] The magnification bias is nonneglible for 
z > 0.5 and large for z > 1 . 



HUBBLE BUBBLE 

The second element of our model consists of using an LTB bubble matched to the EdS background metric of the 
meatball model to describe the local metric around the observer (see Ref. fioll for a discussion about matching in 
swiss-cheese models). The new ingredient is that an observer inside a void expanding faster than the background sees 
an apparent acceleration (see, e.g., Refs. |[Ill[l2tl ). This effect is easy to understand: our cosmological observables are 
confined to the light cone and hence temporal changes can be associated with spatial changes along photon geodesies. 
For example "faster expansion now than before" is simply replaced by "faster expansion here than there". This is why 
a local hubble bubble model can mimic the effect of a cosmological constant at recent times (z < 1). 

The price to pay in a simple hubble bubble model is that the inhomogeneity has to extend to the point in space/time 
where the effect of dark energy disappears. This requires 111 211 an enormous local void of radius 1.5 — 2 /z^'Gpc. 
Moreover, to avoid a too large dipole in the CMB, our position in the void would have to be very special leading to a 
gross violation of the Copernican principle. The model of this letter features inhomogeneities on much smaller scales. 
Here the effects of the dark energy at large redshift (z > 0.5) are mimicked by the lensing bias, and the local void only 
needs to model the cosmological constant at small redshifts where the lensing effects are neglible. 

The spherically symmetric dust LTB metric can be written as 

ds 2 = -c 2 dt 2 + R ' 2{r ;'} 1 dr 2 +R 2 (r,t)da 2 . (9) 
1 — k(r)r z 

This reduces to the usual FLRW metric when R(r,t)/ r = a(r,t) — > a(t) and k(r) — > ±1 or 0. The Einstein equation for 
this metric is 

d 2 (r,t) 8nG c 2 k(r) 

*fit)=— PM-*fri)> (10) 
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FIGURE 5. AH r /H (top panel) and density contrast (bottom panel) with respect to normalized proper distance (R/a(t)) at times 
corresponding to the redshifts indicated in the plot. 



where p is the average density up to the shell r. We specify our LTB model at the recombination time with a density 
contrast of order ~ 10~ 3 in accordance with the CMB. The subsequent evolution of the density contrast is shown in 
Fig. [5] Also shown is the ratio AH r / H where AH r = H r — H, H is the EdS expansion rate and H r = R' / R' is the radial 
expansion rate in the LTB bubble. At late times the initial perturbations generate two fast expanding (AH r > 0) voids of 
radius 50 — 100/z _1 Mpc, surrounded by two thin collapsing (AH r < 0) shells. The inner void of radius ~ lOO/z^'Mpc, 
where the observer is located, has a local expansion rate of h ~ 0.58 which, as anticipated, is larger than the EdS value. 
Moreover, the density contrast of the latter is close to 8 w —0.5 consistently with our meatball model for the overall 
universe. 

In order not to have a too large dipole in the CMB the observer has to be at a distance < 15 hr 'Mpc from the center 
of the inner void. Tighter constrains will come from future probes sensitive to off-center anisotropy (see, for example, 
B). 



RESULTS 

Let us now describe the main features of our results shown in Fig. [6] First, a matter dominated model cannot exhibit 
local acceleration, as one can formally show using the Raychaudhuri's equation |g]. This behaviour is evident in the 
zoom in the bottom panel, where the slope of the inhomogenous model curve is seen to start flat at z = 0. However, 
at finite redshifts the spatial variation of the expansion rate mimics the apparent acceleration and a positive Am 
proportional to the difference between local and background expansion rates Ho — appears. Beyond the dimensions 
of the local hubble bubble, at z > 0.1, the slope of the curve follows the EdS model, until the lensing bias effect becomes 
effective at z > 0.5 and starts to pull the curve up again as is clearly seen in the top panel in fig. [6] 

We performed also a % 2 analysis to estimate the goodness-of-fit of our model. Results are shown in Table [2] 
Obviously the EdS model has a poor fit, but introducing inhomogeneities improves the % 2 significantly. The local 
Hubble bubble has a larger impact on % 2 , but lensing also has a large positive effect. Our model is still worse than 
the ACDM model, but not dramatically so. The main characteristic differences between our inhomogenous model and 
the ACDM model are the smooth rise of the predicted Am at large z in the former model and the different behaviour 
around z ~ 0.1 — 0.2, as seen in Figure [6] Both these issues would be easily resolved by a JDEM-like survey with a 
data set of 2000 SNe. 
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FIGURE 6. Top panel: shown is the distance modulus with respect to the empty universe for the inhomogeneous model of this 
paper, for the ACDM model and for the EdS model together with the full Union Catalogue of Ref. |8]. Bottom panel, zoom for 
small redshifts. 

We would like to stress that we did not explore all the parameter space in order to minimize the x 2 - We indeed chose, 
as explained before, the parameters qualitatively and we calculated the % 2 to give a goodness-of-fit of the meatball 
model with respect to the concordance model. 

Finally let us note that including inhomogeneities changes also the predictions of the ACDM model. In Table|2]we 
give the result of a simulation for the same meatball mass spectrum we used with the EdS model, in the context of a 
ACDM model whose global parameters are set to those used in Ref. [8]. Interestingly, including inhomogeneities to 
the concordance model makes the fit worse and suggests, as discussed in Ref. I14I1 . a smaller value of 

TABLE 2. x 2 for Union Catalogue. 

Model x 2 

ACDM 312 

ACDM + meatballs 323 

EdS 608 

EdS + H. bubble 440 

EdS + H. bubble + meatballs 396 



Of course our setup is but a simple toy model. However, we believe that our results make it explicitly clear that 
different effects of inhomogeneities can pull into the same direction, making the PBS depart from the GBS. Moreover, 
we believe that the two effects of the voids considered here are not the end of the story; to conclude on the viability of 
inhomogeneous universe models as the possible explanation of the apparent acceleration, it is crucial to consider all 
possible effects in as realistic a model as possible. 
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